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Abstract. We give an overview of the phenomenology of strangeness 
enhancement in heavy ion collisions, within the paradigm of the statistical model 
of particle production. We argue that, while strangeness enhancement data is 
suggestive of a phase transition, the mechanism triggering enhancement is still 
elusive. We study the feasibility to constrain this mechanism by determining the 
scaling variable of strangeness enhancement. We further argue that to test the 
applicability of the statistical model generally, and to confirm our interpretation 
of the statistical physics responsible for strangeness enhancement, the scaling of 
strange particle fluctuations (K/tt and other particles) w.r.t. yields has to be 
explored. 
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1. Theoretical motivation 

The objective of heavy ion physics is to produce a locally equilibrated state of quarks 
and gluons, the "quark-gluon plasma" (QGP), and to study the properties of this new 
phase, as well as the transition between quark-gluon plasma and hadronic matter. 

One of the first proposed signatures for the new state of matter was strangeness 
enhancement [TJ [2]. The basic idea was to look at the abundance of strange 
and (especially) multi-strange particles in A-A collisions, and compare with the 
corresponding p — p and p — A collisions (all scaled by the number of participants 
Npart or the multiplicity dN/dn). The appropriate observable is the "enhancement" 
E of strange particle Y 

E (jgj) (^- A ) m . rp (dN/d V ) p „ A (jyg-A) 
Y << Qrt A )(A^) Y (dN/dr,) A _ A ( N r-A^ ^ 

(We shall comment on the two definitions later. Approximately, they give the same 
result, through the corrections could significantly modify the scaling). 

It is expected to that Ey > 1 and increasing with strangeness content of the 
particle. The reason given in p] is that strangeness-producing processes in a QGP 

q q <-> s s g g «-> s s (2) 

should equilibrate faster than the corresponding processes in a hadron gas (HG) 

7T+ n~ <-> K K vr N AK (3) 

this can be seen relatively quickly by computing the momentum exchange (Q) 
needed for these processes ((Q) ~ 2m s for Eq. [2j(Q) ~ tuk for Eq. [3l In 
a thermally equilibrated medium the equilibration time depends on (Q) /T). The 
greater degeneracy of massless quarks and gluons with respect to pions and nuclcons 
makes the difference of thermalization times for the two phases even larger. Hence, 
the strangeness abundance should reach chemical equilibrium (where the strangeness 
relative to light quark abundance depends only on mass difference and temperature, 
not on initial conditions) much faster in a quark gluon plasma than in a hadron 
gas. Since the initial strangeness abundance in collisions is zero, the number of 
strange particles in a system of a certain lifetime with a phase transition should be 
parametrically higher than for a similar system where the transition did not occur. 

This is particularly true for multi- strange particles: In the hadron gas phase, 
these can only be produced sequentially 

TTTrTrN^TrTrAK^TrEKK^ttKKK (4) 

so the equilibration time for these is parametrically -/V s ta. In a hadronizing quark- 
gluon plasma, strange hadrons are presumably produced via coalescence of quarks. 
Hence, A, S and are "automatically" equilibrated to the same degree. 

In summary, strangeness particle abundance in a collision where a QGP is 
produced (a high energy A-A collision) should be enhanced w.r.t. a collision where 
hadronic dynamics is at play (a p-p or p-A collision, or an A-A collision where 
hadronic dynamics dominates). The enhancement should also grow with the hadron's 
strangeness content. The onset of this enhancement could signal the appearance of a 
phase transition, or more generally a change in the degrees of freedom of the system. 

Experimental data have soundly confirmed this prediction [5HH[S]. Both at SPS 
and RHIC, enhancement defined w.r.t. p-p collisions is unity for p-A and d-A. At A-A, 



Phenomenology of Strangeness enhancement in heavy ion collisions 



3 



it increases with both N part and strangeness content of the particle, hiE oc sin N part . 
The constant of proportionality does not vary between RHIC and SPS. This is exactly 
the behavior expected under the hypothesis that, in A- A (but not in p-A) systems the 
density of strangeness is higher than in p-p. 

Given that no evidence of such a scaling violation has been found in other soft 
heavy ion observables (v2 [6j and HBT radii [7] scale remarkably smoothly with 
multiplcitiy) , and the difficulty of doing scaling studies with Hard observables, it 
is worth underlying that the scaling behavior of strangeness enhancement is very 
reminiscent of the behavior associated with a phase transition. If we find an energy 
where enhancement of multi-strange particles does not occur (the coming low-energy 
run at FAIR [9] and the SHINE program at the SPS [TO] will be an optimal 
opportunity to look for such a system) , we would conclusively confirm that strangeness 
is an "order parameter" where scaling violation occurs for a critical energy and system 
size. 

However, the interpretation of this observable has not been universally agreed 
upon. In particular, it was pointed out that the same behavior of enhancement w.r.t. 
strangeness can be obtained through strictly conserving strangeness [TT|. Recently, a 
change in the thermalization time through the excitation of higher hadronic resonances 
and/or multi-particle processes [12] has been suggested to produce the same effects. 

The focus of this work is to review these hypotheses, in particular how strangeness 
enhancement scales with particle type, energy, system size in each. 

2. Statistical hadronization phenomenology of strangeness enhancement 

The idea of modelling the abundance of hadrons using statistical mechanics techniques 
has a long and distinguished history [131 HH HH HBJ. In a sense, any discussion 
of the thermodynamic properties of hadronic matter (e.g. the existence of a phase 
transition) requires that statistical mechanics be applicable to this system ( through 
not necessarily at the freeze-out stage). That such a model can in fact describe 
quantitatively the yield of most particles, including multi-strange ones, has in fact 
been indicated by fits to average particle abundances at AGS, SPS and RHIC energies 
[HI [HI HH [201 [ID H2 [23]. We refer to some of these references for a review of the 
statistical model. For this work, it is sufficient to say that the particle abundance can 
be calculated using the following formula: 

(Ni) = f (V, Q - Q) J] (A^-% Q ^) F (2) + £ b ^ ( Nj ) (5) 

i j—>i 

we now examine each term, explaining its physical origin and its scaling properties. 

The resonance feed-down contribution, Ylj—n bj^i (Nj) (in terms of the Branching 
ratio bj^i and the resonance abundance (Nj)) is straight-forward to take into account, 
provided the full table of resonances is available. In the calculations of this work, we 
have used the data-base in SHARE [22] [23], 

The thermo-statistical weight, in both the Canonical and Grand-Canonical 
ensembles, F (^) , is approximately given in terms modified Bessel function [19] of the 
ratio of the mass of the particle m as well as the chemical freeze-out temperature. Since 
strange particles are heavier than light particles, a higher/lower chemical freeze-out 
temperature could give rise to an effective strangeness enhancement/suppression. It is 
worth underlining that freeze-out temperature is highly correlated, as a fit parameter, 
to the phase space factors 7q examined later [22l [23] . 
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The Volume factor in Eq. [5J / ( V, Q — Q) , is trivially proportional to the volume, 
~ V in the thermodynamic limit (the limit where V/Vq 3> 1, where Vq is the "typical 
scale" of one unit of charge Q) and in the Grand-Canonical ensemble. If strangeness 
abundance is really described in this limit, then all volume effects from the strangeness 
enhancement can be eliminated by defining the enhancement in terms of multiplicity 
(right hand definition of Eq. [IJ, since multiplicity, in the Grand Canonical ensemble, 
tracks the system volume at freeze-out 

Away from the thermodynamic limit (V/Vq < 1), it becomes ensemble-specific. 
In the Canonical and Micro-Canonical ensembles, it will scale, approximately, as the 
volume to the power of the total net charge Q of the particle, F(Q,V) ~ V^^^ 
, since the phase space exactly conserving the charge of a small system producing 
a particle of charge Q is (approximately exponentially with Q) suppressed w.r.t. the 
thermodynamic limit expectation. In the Grand canonical ensemble charge violation 
is allowed event by event, due to the exchange of charge between the system and the 
bath. Thus, the particle yield is proportional to volume for all volumes. 

Because of this effect, if the smaller system is away from the thermodynamic limit, 
the larger system is closer to it and charge is conserved for the observable system [§] 
the Enhancement calculated according to EqQ]is 



Ey ~ / 



N AA 
part 

part 



\Qy-Q y \ 

MAX 



,/(£«!) ocx , lim f (x) — x (Constant) (6) 



This behavior, and the scaling w.r.t. Q — Q are experimental signatures capable of 
ascertaining the contribution of canonical suppression to strangeness enhancement. 

The chemical factor in Eq. O Yl f^Q* < ^ i 7Q* +< ^ i )) counts the contributions of 
the density of the conserved charge inside the system, parameterized by the fugacity 
Aq and the phase space factor jq. The first factor, Aq, counts the net charge density 
of the system, can not give rise to strangeness enhancement (since the net strange 
density is zero independently of strangeness density) and, as it commutes with the 
Hamiltonian, it is independent of the degree of chemical equilibration of the system. 
To account for chemical non-equilibrium, another parameter, jq, can be introduced. 
This counts the abundance of QQ pairs in the system. It is normalized to be = 1 at 
equilibrium and in general, does not commute with the Hamiltonian. 

If strangeness equilibration indeed proceeds faster in the QGP phase, then 7 S is 
expected to be < 1 in p-A systems and ~ 1 in A-A systems. It is then trivial to 
show that enhancement scales as 7f +s . Note that the exponent in the enhancement 
factor is the total strangeness s + s, rather than the net strangeness in the earlier 
canonical suppression scenario (-/ s increases the density of strange particles). For 
most particles, these two numbers are the same, but not for all. For instance, in 4> 
s — s — 0, but s + s — 2. Hence, in the canonical enhancement scenario the <f> should 
not be enhanced at all, but in the scenario where chemical equilibration dominates, 
it should be as enhanced as the H. Other mesons with an ss content, such as the n, 

t Npart tracks the system volume at the start of the evolution. Deviations from isothropic 
expansion, initial state effects, and event-by-event fluctuations of the dynamics 24 a could well make 
the relationship between initial and final volume non-trivial, thereby spoiling any scaling in the 
strangeness enhancement. Because of these ambiguities, we favour defining Ey in terms of the 
multiplicity, through the definition in terms of N par t is more widely used 

§ Not guaranteed at RHIC experiments, where charge could well be exchanged between the observed 
central rapidity region and the "bath" given by the rest of the system. The non-observation of neutral 
particles in most experiments also spoils observable effects of conservation laws 
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should behave accordingly (SHARE 2jJ] takes the mixing of quarks into account in fits 
including 7). Thus, the two scenarios are not different interpretations, but different 
theories, falsifiable by an accurate study of the scaling behavior of enhancement. 

The dependence of strangeness enhancement w.r.t. N part could however be 
similar in both cases: Strangeness enhancement, as we have seen, saturates in the 
large N par t limit since the system there reaches the thermodynamic limit. In the 
"faster equilibration in QGP" scenario, 7 S should also approach the saturating value 1 
asymptotically provided the system obeys a kinetic type evolution equation throughout 
the phase transition and the hadronic phase. In this case, the enhancement is 

\ 7J / N pcLTt — too 

However, there is (possibly) more to this issue: the equilibrium strange particle 
abundance in an ideal QGP should be parametrically higher than the corresponding 
equilibrium abundance of strange quarks in a HG 

3- 1 m m r, Ac 5 >:L m A -m 

~ e T , — ~ e T (8) 

QGP '"1 HG n 1 HG 

Hence, the abundance of strange particles at T > T c should be above that for a 
perfectly equilibrated hadron gas. How does such an abundance evolve when the 
basic degrees of freedom of the system change? Dimensional analysis indicates that 
this can be understood via the dimensionless parameter 



a 



Xs{T, hb) 



Ps(T, p b) 



dT dfii 



dr dr 



(9) 



where x s , Ps are the strangeness susceptibility and equilibrium density in the hot phase, 
7p and are the rates of change in temperature and baryochemical potential in the 
co- moving frame and t s {T,^b) is the equilibration time-scale for s close to the phase 
transition. If a <C 1, the system's chemical composition changes smoothly from the 
old phase to the new phase equilibrium value as the system undergoes the transition. 
If a 3> 1, the system does not have time to adjust its chemical composition, and the 
colder phase is produced out of chemical equilibrium, reflecting the degree of chemical 
equilibrium of the hot phase. In this case, it is expected that 7,, > 1, a value forbidden 
in a Boltzmann equation (and in general unlikely unless a fast change in microscopic 
degrees of freedom occurs) 

For the case of strangeness, Xs/Ps is known from the lattice to have a peak, 
but not a divergence (in accordance to the current understanding of deconfincment 
to be a cross-over, not a first order phase transition). can be estimated through 
hydrodynamics. r s is completely unknown. The chemical equilibration timescale is, 
however, thought to be related to the bulk viscosity, which is presumed to exhibit 
a strong peak [25l |26l [27]- Thus, a might well be considerably larger than a lattice 
estimate, together with a naive extrapolation (e.g. r s ~ f/T 1 ) might suggest. 

The exact same reasoning, however, holds for the light quark degrees of freedom 
q, which could therefore be similarly oversaturated with respect to equilibrium (most 
statistical models assume them to be equilibrated). If this is the case, both j q and 
7 S should be used, and are expected to be > I at freeze-out (their exact value is not 
readily calculable, through it can be estimated by entropy matching [19] ) 

In this scenario, several factors contribute to strangeness enhancement: One is 
greater strangeness content (due to 7^ > I). Another is the number of quarks (7^ > I 
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so Baryons are more enhanced than mesons). The onset of the saturation j q ^ s — » 1 
for a long-lived (N part large) system need not occur 

Putting all these points together, while the interpretation of strangeness 
enhancement is not univocal, different models predict different scaling variables for 
the enhancement. If strangeness enhancement is due to canonical suppression, it 
should increase with net strangeness |s — s\ of the hadron. If it is due to a growth 
of strangeness phase space density, it should grow with total strangeness s + s of the 
hadron. If the phase space density saturates to its equilibrium value up to freezeout, 
strangeness enhancement should also saturate. If, however, lack of equilibrium is due 
to a fast transition from a phase with different equilibrium values for strange and 
light quark density, this saturation does not need to occur. In this latter scenario, the 
strangeness enhancement need not depend entirely on strangeness content, but can 
also show a Baryon/Meson dependence due to light quark non-equilibrium. 

A change with freeze-out temperature with system size (expected if the chemical 
non-equilibrium is associated with supercooling, as in dT5]) could give a further scaling 
of enhancement with particle mass. A non-linear dependence of the final freeze-out 
volume with N part (expected if expansion is not quite isentropic) could lead to an 
enhancement contribution for all particles, independently of mass and chemical content 
(Defining Ey with multiplicity would eliminate this effect). A reinteracting hadron gas 
stage after hadronization could affect resonance enhancement in a non-trivial way [29] 
(through this effect is expected to be ~ FiT^c, where F; is the width of the resonance 
and thg the lifetime of the interacting HG). Since several changes could occur at the 
same time, inverting the observed enhancements to obtain what statistical parameter 
changes between p-A and A-A systems is non-trivial. 

An explicit calculation of enhancements from both the equilibrium and non- 
equilibrium statistical model (parameter fits taken from [28], relevant for RHIC 
energies) is shown in Fig. [T] Summarizing this section, different scenarios for the 
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Figure 1. (Color online) Enhancement factors of different stable particles (left 
panel) and resonances (right panel) for the equilibrium (dashed lines) and non- 
equilibrium (solid lines) statistical models. Red ellipses show where the scaling in 
the models differs 



origin of strangeness enhancement within the statistical model give, generally, different 
scaling observables controlling strangeness enhancement. The numerous competing 
effects, however, render it non-trivial to conclusively establish which model is more 
physically appropriate by scaling studies alone. The applicability of statistical 
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mechanics at all energies and system sizes also remains untested. We now propose 
that the scaling of strange particle fluctuations is necessary to address these points. 

3. Scaling of fluctuations as a confirmation of statistical behavior 

Particle fluctuations are a promising observable to falsify the statistical model and 
to constrain its parameters (choice of ensemble, strangeness/light quark chemical 
equilibrium) [30] . One can immediately see that fluctuations are a stringent statistical 
model test by considering the fluctuation of a ratio between two random variables. 



_ ((AiVi) 2 ) , ((AjV 2 ) 2 ) (ANiAN 2 ) 

°K/N a - {Ni) 2 + {N2) 2 2 {Nl){N2 )- 



Since, for an equilibrated system, in the thermodynamic limit ((AiVi) 2 ) ~ (TV) ~ (V), 
where (V) is the system volume, it is clear that "j^/j^ depends on the hadronization 
volume in a manner opposite to that of particle yields. Volume fluctuations (which 
make a comparison of statistical model calculations to experimental data problematic), 
both resulting from dynamics and from fluctuations in collision geometry, should not 
alter this very basic result since volume cancels out of particle ratios event by event 
|31j . provided hadronization volume is the same for all particles (a basic statistical 
model requirement). 

Thus, observables such as (Ni) Ojvwjv.,) provided (Ni^) and c 2 ^^ are measured 
using the same kinematic cuts, should be strictly independent of multiplicity and 
centrality, as long as the statistical model holds. Furthermore, the value of this 
constant should be calculable with the same parameters used to calculate particle 
yields (T, A 9 . s , 7 g . s and choice of ensemble |33j ) This is a stringent test, because 
variables which generally correlate for yields (such as T and 7 9 , s ) aniz-correlate for 
fluctuations [501 152] . 

The independence of (Aq) o~ 2 Nl f N2 with volume can be stringently tested: if the 
temperature and chemical potentials between two energy regimes are approximately 
the same at freezeout (this should be the case for RHIC top energies and LHC, provided 
chemical equilibrium holds), (N\) cr^ , Na should also stay constant across energies. 
(Ni) ct^i/a^ snou ld also be approximately the same for Au-Au and Cu-Cu collisions 
at high RHIC energy [32], and independent of centrality. If 7 9iS at SPS,RHIC and 
LHC are oversaturated w.r.t. equilibrium (As in the scenario described in [19]), then 
(Ni) 0jVi/iV2 should still De independent of centrality and system size (Cu-Cu vs Au- 
Au) for a given energy range, but should go markedly up for the LHC from RHIC, 
because of the increase in 7 9 and 7 S . 

While fluctuations are more sensitive to acceptance cuts then yields, these 
systematic errors can be at least partially removed via mixed event subtraction 
[34] . An appropriate observable to model and measure [34] [35] [36] is therefore 
a dyn = cr 2 — °~mix where is the mixed event widtl{jj]. In the absence of any 

correlations, a mix reduces to the Poisson expectation, o-^ nix = 1/ (Ni) + 1/ {N 2 ). We 
therefore propose to use ^j^i/iVa = a dyl N2 to tes t the statistical model among 
different energy, system size and centralities. 

We have calculated ^%\/ff 3 f° r RHIC and LHC energies, for the sets of parameters 
used in [37]. The left and right panel in Fig. [5] shows what effect three different sets 

| the Udyn observable, equivalent to o\ for large systems [34] . is more appropriate for rarer particles 
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Figure 2. Scaling of the observable w.r.t. multiplicity for particle ratios at 
RHIC (left panel) and LHC (right panel) energies in the ■y qyS > 1 non-equilibrium 
scenario, equilibrium scenario and equilibrium-canonical scenario 



of j q , s inferred in [37 would have on ^^-i^- and ^\-/ K +- In the left panel we 
have also included the value of i^- for top energy RHIC. As shown in [32] . this 
value for top centrality matches expectations for non-equilibrium freeze-out (and is 
significantly above equilibrium freeze-out). A centrality scan of i v - , necessary to 
confirm the consistency of this result has not, however, as yet been performed. 

The scaling should break if global correlations persist, such as is the case if 
the Canonical and micro-canonical ensembles |33j are physically more appropriate 
to describe the system (and Canonical suppression contributes to strangeness 
enhancement). If strangeness at RHIC/the LHC is created and maintained locally, 
^n 1 1 /n 2 snou ld develop a "wiggle" at low centrality, and be lower (by half for ^^+/ K - ) 
than Grand Canonical expectation. 

In conclusion, we have outlined the different statistical models capable of 
explaining the experimentally observed onset of strangeness enhancement, and argued 
that the scaling of enhancement with hadron type is generally different in each of these 
models. We have further argued that all statistical models can be stringently tested 
by observing the scaling of strange particle fluctuations w.r.t. yields. 
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